Boer-Mulders function of the pion in the MIT bag model 
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We apply the MIT bag model to study the Boer-Mulders function of the pion, a T-odd function 
that describes the transverse polarization distribution of the quark inside the pion. We simulate 
the effect of the gauge link through the "one-gluon-exchange" approximation. We consider both the 
quark helicity nonflip and double-flip contributions. The result in the MIT bag model is compared 
with those in the spectator models. 
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One of the main tasks in QCD and hadron physics is to 
understand the transverse partonic structure of hadrons, 
especially the nucleon and the pion. The inclusion of 
the parton transverse motion introduces new types of 
parton structure, the so-called transverse momentum de- 
pendent (TMD) distributions, or alternatively the three- 
dimensional parton distribution functions in momentum 
space. They extend the concept of traditional Feynman 
distribution functions and encode a wealth of new infor- 
mation on the nucleon structures [2-0] that cannot be de- 
scribed merely by the leading-twist collinear picture. Of 
particular interests are the leading twist T-odd TMD dis- 
tribution functions, such as the Sivers function [Ullj and 
the Boer-Mulders function [8j. They arise from the cor- 
relation between the nucleon/quark transverse spin and 
the quark transverse momentum, and they can account 
for the polarized and unpolarized spin asymmetries in the 
semi-inclusive deeply inelastic scattering (SIDIS) 
and the Drell-Yan |15l - [l7j processes. 

As a spin-0 hadron, the pion has a simpler par- 
tonic structure than that of the nucleon, i.e., in lead- 
ing twist there are two TMDs of the pion, the unpo- 
larized TMD fi(x, fe T ) and the Boer-Mulders function 
hi(x, fe|,). However, the pion TMDs are less known 
from experiments than those of the proton since they 
cannot be probed in the SIDIS. Theoretically, the Boer- 
Mulders function of the pion has been studied by lattice 
calculation [l8| and model calculations [l9r]24j|. In the 
latter case, different treatments on the gauge link have 
been used, namely, the one-gluon exchange approxima- 
tion JSEl [24| and the nonperturbative eikonal meth- 
ods |22l.l23j. which take into account higher order gluonic 
contributions, respectively. In this paper, we study the 



Boer-Mulders function of the pion using an alternative 
model, the MIT bag model [25(. This model has been 
applied to study the TMDs of the proton, including the 
T-even distributions (26|, the Sivers functions [U{2^], 
and the Boer-Mulders functions [13, The calcula- 

tion of the T-odd TMDs by the MIT bag model has pro- 
duced their main features, for instance, the sign and the 
Burkardt sum rule [3l[ for the Sivers function (29|, and 
the sign for the Boer-Mulders function [32j]. Therefore, it 
is worthwhile to use the same model to study the Boer- 
Mulders function of the pion. 

Unlike the Boer-Mulders function of the proton, which 
can be probed in both the SIDIS process and the Drell- 
Yan process, the Boer-Mulders function of the pion may 
only be detected in the Drell-Yan process. Fortunately, 
the new ttN Drell-Yan program will be conducted by 
COMPASS HI at CERN very soon; also there is a irN 
Drell-Yan plan proposed by SPASCHARM [||. The up- 
coming Drell-Yan experiments can achieve unpolarized 
and polarized scattering, so they will provide the oppor- 
tunities [35|, [36( to access the chiral-odd TMDs of the 
pion as well as the nucleon. 

The quark-quark correlation function for the pion has 
the form 
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is the gauge link (Wilson line) connecting the two dif- 
ferent space-time points and £ by all possible ordered 
paths followed by the gluon field A running along a 
process-dependent path. In this work we calculate the 
Boer-Mulders function in the SIDIS process. 
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The leading-twist TMDs of the pion can be obtained The TMD distribution f%(x, kj,) can be calculated 
from the correlator fcy) by the following traces: straightforward in the MIT bag model, in which the 

quark fields are expressed in the following general form: 



e ap k 



f 1 (x,k 2 T ) = ±Tr[<S>(x,k T ) 1 +] 1 (3) 
M ^hi{x, k 2 T ) = ^Tr k T )ia +a l5 ] . (4) 



i S a (x,t) = ^ N{b a {nKm)ipnKjm(x,t) + d^(nKm)ip- n - K j m (S,t)} (5) 

n>0,K=±l,m=±l/2 



where ip is the wavefunction in the position space. After 
performing the Fourier transformation, one obtains the 
momentum space wavefunction of the quark (25j . 



to(k) Xr 



tp m (k) = WAnNRl 



(J ■ fcti(fc) Xr, 

and the wavefunction of the antiquark, 
<f>m(k) = 



(6) 
(7) 

where k = j is a unit vector with k = |fe|, oj w 2.04 for 
the lowest mode, Rq the bag radius, Xm the Pauli spinor, 
a the Pauli matrix, and N the normalization factor with 
the form 



N 



2R 3 Q (w - l)j 2 {u) 
The functions U(k) are calculated from 
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where j% are the spherical Bessel functions. 

Using the isospin symmetry and charge-conjugation 
operation, the unpolarized TMDs of the charged pion 
can be connected by 



h = h = h — h = /It- l iU J 

The function /i w can be calculated by inserting the quark 
field in the MIT bag model into the correlator ([T]) in the 
absence of the gauge link 



hw(x, k T ) 
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where k z = %, and k z = xM % — e with e — uj/Rq. The 
distribution for the neutral pion is a half of . 



The MIT bag model has also been extended to calcu- 
late T-odd TMDs [H HI, [H , such as the Sivers function 



k, m.2 








q 







w ^ • 

k 3 - q, m 4 


&3 , TU3 



k — q, mi 



+h.c. 




FIG. 1: One gluon exchange contribution to the pion 
Boer-Mulders function. The graph has been drawn using 
JAXODRAWlH. 



and the Boer- Mulders function of the nucleon. As in the 
original MIT bag model there is no explicit gluon degree 
of freedom, which is crucial for nonzero T-odd TMDs; in 
these calculations the effect of the gauge link is incorpo- 
rated by introducing "one-gluon-exchange" [13, [111 H3] 
or invoking instanton effects (28|. In our calculation of 
the pion Boer-Mulders function, we follow the former ap- 
proach to expand the gauge link to order O(g) to obtain 
the expression 
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where we have used the covariant gauge. We also point 
out that the calculated distribution is for the semi- 
inclusive DIS process. The corresponding diagram is 
shown in Fig. [TJ In Eq. (IT2l we use T a and T a (the con- 
jugate representation) to denote the Gell-Mann matri- 
ces associated with the quark and the antiquark, respec- 
tively. They are related by 



-(T°) 



a\T 



(13) 



Defining 

Cm 1 m 2 ,m. 3 m 4 , = 2y ^jy {P^ | b mi &^j 2 rf TO3 d mi \Ptt) j (14) 

we obtain the following nonzero spin coefficients: 

2 2 

C+ + - = C ■++ = C + _,_ + = C- +>+ 
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These coefficients have already included the contribution 
from the color factor 



(16) 



ijkl 



calculated from the color (-singlet) structure of the pion 
and the Gell-Mann matrices. 

Inserting the bag wavefunctions © and (JTJ), and the 
spin coefficients (|15l) into (|T2|) . we arrive at the final ex- 
pression of the pion Boer-Mulders function 



hi,(x,k 2 T ) = -g 2 s ^- J |^rj Jd 3 k 3 (R ++ R' ++ + I ++ I' ++ +R + -R' + _+I + -I' + _), 



(17) 



where 
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In Eq. (JT7|) the real functions R mi m 2 , I mi m 2 , R' m „ m ,, and 



^m 3 m 4 are defined as 
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There are totally 16 functions, among which 8 are in- 
dependent and have the following forms: 
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R++ = kJoik'Mk) - k' x t (k)h{k') + (k' z k - ( 21 ) 

R+- = -t (k')t (k) - k z t (k%(k) - k'Mk)ti{k') + {-k% - k' x k x + kykyMk'^k), (22) 

I++ = k y t (k')h(k) + k' y t (k)h(k') + (k' z ky + k'ykMk'Mk), (23) 

I+- = CKh + k' y k)h(k')h(k), (24) 

R++ = k ■ k' 3 ti{k3)ti(k' 3 ) + kzti(k 3 )t a (k' 3 ) +k' 3z t (k3)ti(k' 3 ) + t (k 3 )t (k' 3 ), (25) 

R+- = -(kJ' 3z ~ kzkxMksMk) - kMk 3 )t Q {k' 3 ) + k' 3 xto{k 3 )h{k' 3 ), (26) 

I'++ = (k 3x k' 3y ~ kyk' 3 Mk 3 )ti(k' 3 ), (27) 

I' + _ = (k 3 yk' 3z ~ kzkyMkMk'z) + kytl{h)t {k' 3 ) - fc^ofe)*!^), (28) 



where k' = \k - q T \, k' = ^ 2L , and k' 3 = |fc 3 - q T \, 
k' 3 = -^t 221 - The functions listed in Eqs.flH])-® agree 
with the functions ii ; 2, Fi,2, Hi,2, and Ji,2 listed in the 
appendix of Ref. |30| . The first two terms on the right- 
hand side of (|17l) are the quark helicity nonflip contribu- 
tions, while the last two terms are the contributions re- 
ceived from the helicity double-flip of quarks. An impor- 
tant observation in the bag model calculation of the pro- 
ton Sivers function [29[ and Boer-Mulders function [30| 
is that, apart from the helicity nonflip contributions, the 
double- flip terms (especially the R'^ term ) is sig- 
nificant and should not be ignored. In the light of this 
finding, in this work we consider both the helicity nonflip 
and double-flip contributions to the pion Boer-Mulders 
function. 

To give a numerical estimate of the pion Boer-Mulders 
function in the MIT bag model, we need to fix the pa- 
rameters in the model, especially the bag radius Rq of 
the pion. In the calculation of the proton TMDs (27l . [29l . 
I30L Wl\ the bag radius is determined by the relation |38l ] 

„ 4 nuj 

*° = 3^' (29) 

where n is the quark (antiquark) number in the bag. 
Here we use the same ansatz for the bag radius of the 
meson. For the str ong coupling a s , we follow the choice 
a s /(47r) = 0.13 in [30|, where the same model has been 
used to calculate the proton Boer-Mulders function. To 
get the appropriate tendency of the distribution at the 
region x — > 1, we use the constraint 5(1 — x — x 3 ) when 
performing the integration in (|12[) , where x 3 = k 3 / P + . 

The left panel of Fig. [5] shows the first fcf, moment of 
the Boer-Mulders function hf^\x), which is defined as 

hfj 1 \x) = |d 2 fe T (A)^(x,^). (30) 

The solid and the dashed curves represent the total re- 
sult and the result contributed by the quark helicity non- 
flip terms. The comparison of these two curves indicates 
that the helicity nonflip and double-flip contributions are 



equally important to the pion Boer-Mulders function. 
Our results show that hj^ is negative, in agreement with 
spectator model and lattice calculations. The sign of the 
pion Boer-Mulders function is also consistent [32| with 
the sign of the Boer-Mulders functions of the nucleon in 
the MIT bag model. 

In the right panel of Fig. [5] we compare xhf£ (x) in 
the MIT bag model (shown by the solid line) with that 
in the spectator model (shown by the dashed line) [l9j . 
where the one-gluon exchange approximation is also used. 
When obtaining the two curves in the right panel of Fig. [5] 
we use the same strong coupling a s /(47r) =0.13 for com- 
parison. The size of xhf^\x) in the MIT bag model is 
smaller than that in the spectator model, while the x de- 
pendence of the distribution are similar in both models; 
that is, they peak at the region x ~ 0.5. 

We also point out that the size of xh^ X \x) in our 
calculation is comparable with the spectator model cal- 
culation that employs the nonperturbative eikonal meth- 
ods [23|. The x dependence of xh^ 1 \x) in these 
two different calculations differ from each other since in 
[111, [23| the distribution peaks at x ~ 0.2. 

Similar to the MIT bag model calculations for the pro- 
ton TMDs, our calculations are performed at the low en- 
ergy bag scale, and we have not considered the evolution 
effect during the entire calculation. Recently substantial 
progress [33, H3| on the evolution of the proton TMDs 
has been achieved. A potential issue for future study is 
to investigate if the same approach can be applied to the 
pion TMDs, especially the Boer-Mulders function, where 
the MIT bag model calculations could be the initial in- 
puts at the low energy scale. Only with full knowledge 
of the initial inputs and the evolution of the pion TMDs, 
we can get more precise predictions of the experiments. 

In summary, we have applied the MIT bag model to 
study the TMDs of the pion. Particularly, we calculated 
the pion Boer-Mulders function, which is a T-odd chiral- 
odd distribution. To obtain a nonzero result, the effect 
of the gauge link is simulated by introducing the "one- 
gluon-exchange" effect. We consider both the helicity 
nonflip and double-flip contributions to the pion Boer- 
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FIG. 2: Left panel: The first fcy moment of the pion Boer-Mulders function in MIT bag model. The solid line corresponds 
to the full result including both the helicity nonflip and double-flip contributions. The dashed line corresponds to the result 
with only the helicity nonflip contribution. Right panel: Comparison of xh^~ (x) in the MIT bag model (solid line) and the 
spectator model (dashed line). 



Mulders function. We estimated the pion Boer-Mulders 
function numerically, showing that it is negative in the 
MIT bag model. We compare our result with the avail- 
able spectator model calculations. Our study provides 
further knowledge on the transverse parton structure of 
the pion. 
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